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The transition probabilities describing the evolution of a neutrino with a given energy along some 
ray through a turbulent supernova are random variates unique to each ray. If the source of the 
neutrinos were a point then all neutrinos of a given energy and emitted at the same time which 
were detected in some far off location would have seen the same turbulent profile therefore their 
transition probabilities would be exactly correlated and would not form a representative sample of 
the underlying parent transition probability distributions. But if the source has a finite size then 
the profiles seen by neutrinos emitted from different points at the source will have seen different 
turbulence and the correlation of the transition probabilities will be reduced. In this paper we 
study the correlation of the neutrino transition probabilities through turbulent supernova profiles 
as a function of the separation Sx between the emission points using an isotropic and an anisotropic 
power spectrum for the random field used to model the turbulence. We find that if we use an 
isotropic power spectrum for the random field, the correlation of the high (H) density resonance 
mixing channel transition probability is significant, greater than 0.5, for emission separations of 
5x = 10 km, typical of proto neutron star radii, only when the turbulence amplitude is less than 
C* ^ 10%; at larger amplitudes the correlation in this channel drops close to zero for this same 
separation of Sx = 10 km. In contrast, there is significant correlation in the low (L) density resonant 
and non-resonant channels even for turbulence amplitudes as high as 50%. Switching to anisotropic 
spectra requires the introduction of an 'isotropy' parameter ki whose inverse defines the scale below 
which the field is isotropic. We find the correlation of all transition probabilities, especially the H 
resonance channel, strongly depends upon the choice of ki relative to the long wavelength radial 
cutoff ki, . The spectral features in the H resonance mixing channel of the next Galactic supernova 
neutrino burst may be strongly obscured by large amplitude turbulence when it enters the signal 
due to the finite size of the source while the presence of features in the L and non resonant mixing 
channels may persist, the exact amount depending upon the degree of anisotropy of the turbulence. 

PACS numbers: 47.27.-i,14.60.Pq,97.60.Bw 



I. INTRODUCTION 



The neutrino signal from the next core-collapse super- 
nova in our Galaxy will give us an unprecedented op- 
portunity to peer into the heart of an exploding star 
and confront our current paradigm of how these stars 
explode with observations. But decoding the message 
will be no easy task because the neutrino signal will have 
experienced so many flavor-changing events on the trip 
from proto-neutron star to our detectors that scramble 
the information, see for example Kneller, McLaughlin & 
Brockman [Tj] and Lund & Kneller 0]. The first flavor 
changing effect the signal experiences is due to neutrino 
self interactions / collective effects in the region up to 
~ 1000 km above the proto neutron star [3l-[l9| followed 
by the Mikheyev, Smirnov & Wolfenstein (MSW) [Him 
effect which is complicated in supernovae by the impres- 
sion of the Shockwave racing through the stellar mantle 
li, 2, 22 27]. Turbulence in the mantle, seeded during the 
earlier neutrino heating/Standing Accretion Shock Insta- 
bility phase [28l - [35j , also needs to be included usually by 
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modelling [1^ l36M39f . Finally, there is the possibility 
of Earth matter effects leaving an imprint in the signal 
though a recent study expects this effect to be minimal 

What makes decoding the signal even more of a chal- 
lenge is that the neutrinos we receive at a given instant 
and with a given energy will not have experienced the 
same flavor evolutionary history. The neutrinos arriving 
at a detector will have been emitted from different loca- 
tions at the source and both the neutrino collective and 
the MSW-|-turbulence effects will vary from trajectory to 
trajectory. Starting with Duan et al. [5|] , the self inter- 
action effects in calculations where the neutrino emission 
over the source is assumed to be spherically symmetric 
have been seen to be 'angle dependent' in the sense that 
a neutrino following a pure radial trajectory differs from 
one emitted at an angle relative to the normal. Presum- 
ably allowing for aspherical source emission would only 
make the trajectory dependence even stronger. Similarly 
the MSW plus turbulence effects are also trajectory de- 
pendent. If we temporarily cast aside the turbulence and 
focus on the gross structure of the explosion i.e. the low- 
est angular multipole moments, an aspherical passage of 
the shock through the star, by itself, leads to a line-of- 
sight dependence. But, one must recall that we will not 
observe the neutrinos from a supernova at widely differ- 
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ent lines of sight, all our detectors are here on Earth. The 
size of the source is of order the proto-neutron star ra- 
dius, i.e. ~ 10 km while the shock effects show up in the 
signal when the shock has propagated out to r ^ 10"* km. 
As long as the curvature of the shock is over a lengthscale 
greater than the source size the neutrinos which appear in 
our detectors will all have seen essentially the same pro- 
file. When one re-inserts the turbulence into the profile, 
one realizes this approximation may no longer be valid 
because turbulence extends to much smaller lengthscales 
even when the shock is far out in the stellar mantle. The 
density profile along two, essentially parallel lines of sight 
to a distant detector separated by ^ 10 km will no longer 
be negligibly dissimilar and one must consider how the 
dissimilarity of the profiles propagates to the neutrinos. 
Any correlation will lead to a potential new feature of 
the neutrino signal. 

It has been shown that the transition probabilities for 
a single neutrino - the set of probabilities that relates 
the initial state to the state after passing through the 
supernova - is not unique when turbulence is inserted 
into a profile: it will depend upon the exact turbulence 
pattern seen by the neutrino as it travelled through the 
supernova [s^ l39l , |4l| . Those transition probabilities are 
drawn from distributions whose properties will depend 
upon the stage of the explosion, the character of the tur- 
bulence, and the neutrino energy and mixing parameters. 
If the coherence of two neutrinos emitted at the same 
time and with the same energy but from different loca- 
tions is small then the final states are uncorrelated and 
one would expect that the flux at a detector would just 
be the mean of whatever distribution describes the tran- 
sition probabilities multiplied by the initial spectra. But 
if the coherence is high then all the neutrinos will have 
the same set of transition probabilities which one might 
expect to 'scintillate' together as the turbulence evolves. 
Of course, this ignores the issue of energy resolution and 
temporal binning of the signal that becomes necessary 
because of the limited statistics. 

The purpose of this paper is to consider the issue of fi- 
nite source size and the correlation of the neutrino transi- 
tion probabilities along parallel trajectories through tur- 
bulent supernova profiles. Our calculations expand upon 
the work of Kneller & Volpe [11] and Kneller & Mauney 
[ilj upon which we rely heavily for the techniques used 
to calculate the turbulence effects and as context for our 
results. We begin by describing the calculations we un- 
dertook paying particular attention to the construction 
of the random fields used to model the turbulence. The 
basic approach to determining the effects of turbulence 
are then demonstrated, followed by the computation of 
the transition probability correlation as a function of the 
separation between the emission points. We finish by 
summarizing our findings and discuss the implications 
for the Galactic neutrino burst signal. 



II. DESCRIPTION OF THE CALCULATIONS 

The neutrino transition probabilities are the set of 
probabilities of measuring some neutrino state Vi given an 
initial neutrino state Vj \.e P{vj ^ Vi) = Pij. We shall 
denote antineutrino transition probabilities by Pij. If 
the S'-matrix relating the initial and final wavefunctions 
is known then these probabilities are just the square am- 
plitudes of the elements of S. The 5-matrix is calculated 
from the Schrodinger equation 



where H is the Hamiltonian. The Hamiltonian is the sum 
of the vacuum contribution Hq and the MSW potential V 
which describes the effect of matter. The vacuum Hamil- 
tonian is diagonal in what is known as the 'mass' basis 
and in this basis Hq is defined by two mass squared dif- 
ferences Smfj = — rrij and the neutrino energy E. 
The mass basis is related to the flavor basis by the Maki- 
Nakagawa-Sakata-Pontecorvo [13, unitary matrix U. 
The most common parametrization of U is in terms of 
three mixing angles, 9i2, Ois and O23, a CP phase and 
two Majoranna phases. 

The MSW potential V is diagonal in the flavor basis 
because matter interacts with neutrinos based on their 
flavor. The neutral current interaction leads to a contri- 
bution to V which is common to all flavors. This may 
be omitted because it leads only to a global phase which 
is unobservable. The charged current potential only af- 
fects the electron flavor neutrino and antineutrinos and 
is given by \/2G'F"-e(r) where Gf is the Fermi constant 
and ne(r) the electron density. 

In matter the two contributions to H means neither 
the mass nor the flavor states diagonalize the matrix. But 
there is a basis known as the matter basis which does di- 
agonalize H i.e. for a given value of the electron density 
there is a matrix U such that WHU = K where K is the 
diagonal matrix of eigenvalues. When the MSW poten- 
tial vanishes the matter basis becomes the mass basis up 
to arbitrary phases. The matter basis is the most useful 
for studying the evolution of neutrinos through matter 
because it removes the trivial adiabatic MSW transition 
and it will be the basis we use to report our results in 
this paper. We refer the reader to Kneller & McLaughlin 
(4^ and Calais, Kneller & Volpe [13] for a more detailed 
description of the matter basis. 

We now turn our attention to the turbulent density 
proflles through which we shall send our neutrinos. As 
usual, we shall model the turbulence by multiplying a tur- 
bulence free density profile by a Gaussian random field. 
Since the spatial extent of the neutrino emission, of or- 
der 10 km, is much smaller than the radial location of the 
turbulence, of order r ~ 10^ — 10^ km, we shall ignore 
any curvature of the density profile features and use a 
plane-parallel model for the supernova. The z axis of our 
Cartesian co-ordinate system is aligned with the radial 
direction of the profile. The profile we adopt is from a 
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FIG. 1: The turbulence free MSW potential as a function 
of distance through a supernova taken from a hydrodynam- 
ical simulation. The vertical lines indicate the positions of 
the reverse and froward shock in the profile. The horizontal 
dashed-dotted line is the two-flavor resonance density for a 
25 MeV neutrino with mixing angle sin^ 26 = 0.1 and mass 
splitting 5m^ = 3 x 10"^ eV^ 



one-dimensional hydrodynamical simulation of a super- 
nova taken from Kneller, McLaughlin & Brockman 
This profile is shown i n fig ure ([T]) and is the same one used 
in Kneller & Mauney [41| . The figure shows the presence 
of two shocks: the forward shock at Ts and the reverse 
shock at rr- In multi-dimensional simulations of super- 
nova both these shock fronts are aspherical and fluid flow 
through the distorted shocks leads to strong turbulence 
in the region between them. Our selection of this profile 
also determines the neutrino energy we shall use since we 
wish the neutrinos to have an H resonance density that 
does not intersect the shocks. Therefore we pick 25 MeV 
for the neutrino energy and the reader may observe that 
the two-flavor resonance density for a 25 MeV, shown in 
the flgure, does not intersect the shocks as required. 

The turbulence is inserted by multiplying the profile in 
the region between the reverse and forward shocks by a 
factor 1 + F{r) where F{r) is a three-dimensional Gaus- 
sian random field with zero mean. The random field is 
represented by a Fourier series, that is 



F(r) = a tanh ( ] tanh 



{An COS (k„ ■r)+ Bn sin (k„ • r)} . 



(2) 



In this equation the parameter sets the amplitude 
of the fiuctuations while the two tanh functions are in- 
cluded to suppress fiuctuations close to the shocks and 
prevent discontinuities. The parameter A is a damping 
scale which we set to A = 100 km. The random part of F 
appears in the the second half of equation ([2]) because the 
set of co-efficients {A} and {B} are independent standard 
Gaussian random variates with zero mean. The k„ are 
a set of wavenumbers and, finally, the paramaters Vn are 
k-space volume co-efficients. The method of fixing the 
Nk fc's, F's, A's and B's for a realization of F is 'variant 



C' of the Randomization Method described in Kramer, 
Kurbanmuradov, & Sabelfeld |45|] which we have general- 
ized to three dimensions. The Randomization Method in 
general partitions the space of wavenumbers into Nk re- 
gions and from each we select a random wavevector using 
the power-spectrum, E(k), as a probability distribution. 
The volume paramaters Vn are the integrals of the power 
spectrum over each partition if the power spectrum is 
normalized to unity. Variant C of the Randomization 
Method divides the k-space so that the number of par- 
titions per decade is uniform over A^^ decades starting 
from a cutoff scale ki^. Throughtout this paper we shall 
use a wavenumber cutoff k^, set to twice the distance be- 
tween the shocks i.e. k^, — n/{rs ~ ?>). The logarithmic 
distribution of the modes increases the efficiency of the 
algorithm in the sense that we can use a 'small' value 
of Nk and also the agreement between the exact statis- 
tical behavior of the field and that of an ensemble of 
realizations is uniform over some range of lengthscales 
i.e. it is scale invariant. This feature is important for 
our study because the oscillation wavelength of the neu- 
trinos is constantly changing as the density evolves. The 
minimum lengthscale we need to cover has been shown 
by Friedland & Gruzinov [s^ and Kneller, McLaughlin 
& Patton [i^ to be the reduced oscillation wavelengths 
for the neutrinos and antineutrinos i.e. A^ = l/|<5A:y | 
and Xij — l/\6kij\ - where Skij and Skij are the differ- 
ences between the eigenvalues i and j of the neutrinos 
and antineutrinos respectively. Kneller & Mauney [4]| 
showed the wavelengths in the turbulence region were of 
order 1 km or greater which is approximately four orders 
of magnitude smaller than the shock separation. This 
means we need to pick A^^ > 4 to cover the necessary 
decades in k-space. 



A. The power spectrum 

The final component of our calculations we have yet 
to discuss is the power spectrum E(k.). In this paper 
we shall consider two power spectra and our first choice, 
due to its simplicity, is a normalized three-dimensional, 
isotropic inverse power-law spectrum given by 



(a - 1) / fc* 



e{\k\-k,). (3) 



for |fc| > k-i, where |fc| is the magnitude of the wavevector 
k. Throughout this paper we shall adopt the Kolmogorov 
spectrum where a = 5/3. The one dimensional power 
spectrum for the fc^ component of the wavevector is 
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/ k^ Y 



2aki, \ |fc; 
2afc+ 



Q{\kz\-K) 

e(fc.-|fc,|) (4) 



which differs from the one dimensional power spectrum 
used by Kneller & Mauney ^ij because for \kz\ > fc^ 
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the power is suppressed by the factor l/a and the one 
dimensional spectrum is non-zero for \kz\ < fc^. The two- 
point correlation function B{5r) for this choice of a power 



where T{n,x) is the incomplete Gamma function. There 
is one last quantity to determine: the number of Nk of 
elements in the sets of random wavenumbers, coefficients 
and volumes. To find this quantity we compare the sta- 
tistical properties of an ensemble of random field realiza- 
tions with the exact expressions as a function of the ratio 
Nk/Nd for a given Nd- The statistical property we com- 
pute is the second order structure function G2 ((5r) which 
is given by 

G2{5y) ^ {F{y + 5y) - F{Y)f (6) 

where 5y is the separation between two points. The func- 
tion G2{5y) is related to the two-point correlation func- 
tion B{Sy) via G2{Sy)/2 = 1 — B{6y). For the isotropic 
power spectrum both G2 and B are only functions of 
the magnitude of Sy and the correlation function is given 
above. In figure ([2|) we show the ratio R{Sr) of the numer- 
ically calculated structure function for the isotropic ran- 
dom field to the exact solution as a function of the scale 
ki,Sr when we use either Nk = 50 wavenumbers spread 
over Nd = 5 decades or Nk = 90 wavenumbers over 
Nd — 9 decades. The numerical calculation is the aver- 
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FIG. 2: The ratio of the structure function G2(5r) as a func- 
tion of fc*5r for two randomly orientated points in a 3-D ho- 
mogeneous and isotropic Gaussian random field to the exact 
structure function. The two curves in the figure correspond 
to {Nk,Nd} = {50,5} (blue solid) and {Nk,N4 = {90,9}. 
At every fc*5r we generated 30, 000 realization of the field and 
the error bar on each point is the standard deviation of the 
mean F{r + Sr) — F{r).. 



spectrum depends only the magnitude of the separation, 
(5r, and may be calculated analytically to be 



(5) 

I 

age of 30, 000 realizations of the turbulence and the error 
bar on each point is the error on the sample mean. The 
figure indicates that the method we use to generate ran- 
dom field realizations reproduces the analytic results for 
the structure function very well and with high efficiency 
because good agreement between the statistics of the en- 
semble and the exact result requires just Nk/Nd = 10. In 
fact we find even Nk/Nd ratios of just Nk/Nd ~ 2 — 3 are 
sufficient to give acceptable agreement but we re-assure 
the reader we shall stick with Nk/Nd = 10. 

But isotropic and homogeneous three-dimensional tur- 
bulence is perhaps not a realistic scenario for supernova 
because the gravitational potential and the general fluid 
flow are in the radial direction. Only on sufficiently small 
scales should the turbulence become isotropic. This divi- 
sion into large and small lengthscales indicates wc should 
partion the power-spectrum so that for | fc^ | > fc/ the 
spectrum is isotropic, where fc/ is the isotropy scale, be- 
tween ki, < \kz\ < ki the spectrum is anisotropic and 
then below the cutoff scale, \kz\ < A;^,, the power spec- 
trum should be set to zero since there should be no modes 
on scales larger than 1/fc*. For \kz\ > kj where the spec- 
trum is isotropic we use a power spectrum resembling 
equation ^ 

Note the additional factor of a in the numerator. For 
ki, < \kz\ < kj we write the spectrum as the product 
E{kx, ky,kz) = E{kx,ky)xE{kz)- The spectrum i^(fcz) is 
chosen to be a continuation of the inverse power-law given 
above while the spectrum in the xy directions, E(kx, ky), 
is the spectrum of the isotropic/homogeneous region in 
these directions fixed at \kz\ = fc/. The spectrum for 
k* < \kz\ < kj is thus 

E(k k k)- f y ( x"^' 

Q{ki~\k\)Q{\k\-K). (8) 

The reader may verify the power spectrum defined 
by equations (O and ([5]) is normalized. This 
anisotropic three-dimensional power spectrum yields a 
one-dimensional spectrum along the z direction given by 

i^i(M = ^(^)"e(|..|-M (9) 
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B{5r) = ' ~ (27rA:^^r)""^ |exp(^) r(-a, 2z7r /c^ Jr) - exp (^) r(-a, -2z7r fc^ (5r)} 
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for \kz\ > ki, which is exactly the same as the one- 
dimensional spectrum used in Kneller & Mauney (4]| . 
There is no analytic formula for the two-point structure 
function for randomly orientated separations using this 



power spectrum but if we consider the two-point struc- 
ture function of the random field for points orientated 
along the z direction then we can compute that in this 
direction 



B{5z) ^ ^ (27rfc^(5z)° ^ |exp r(l - a, 2z7r fc* (5z) -|- exp (^^) T[l - a,-2nr Sz)^. (10) 



Compared to the isotropic spectrum above, this 
anisotropic spectrum differs in important ways. First, 
even if we set kj — fc^ we observe that the lack of power 
in the region \kz\ < fc^ means we have to compensate 
by increasing the structure / decreasing the correlation 
by the factor a. This increase is the reason for the ap- 
pearance of the extra factor a in equation ([7]). Next, 
as we increase the ratio // = kj/ki,, we push more and 
more of the structure of the field in the xy direction to 
ever smaller scales reducing even further the correlation 
of the field at some fixed non-radial separation 5x com- 
pared to the isotropic case. This extra power at small 
scales can be seen in figure ([3]) which is a plot of the ra- 
tio of the one-dimensional two-point structure function 
in the x direction relative to the structure function along 
the z direction at the same separation scale for three val- 
ues of //. As promised, when // = 1 there is an equal 
amount of structure in the field along both radial and 
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FIG. 3: The ratio of the structure function G2{Sx) as a func- 
tion of ki,5x for two points aligned along the x direction to 
the structure function G2{Sz) of two points aligned along the 
z direction at Sx = Sz. The three curves in the figure corre- 
spond to ki = fc* (solid), ki = 10 fc* (dashed) and kj = 100 fc* 
(dot-dashed). At every k^,Sx we generated 30,000 realization 
of the field and the error bar on each point is the standard 
deviation of F{x + Sx) — F{x). The structure function G2{Sz) 
was computed using the correlation function given in (|10p 
and the relationship G2{Sz)/2 = 1 — B{Sz). The inputs to 
the random field generator were Nk ~ 90, Nd = 9. 



non-radial directions but as // increases we push more 
and more of the structure of the field in the xy direction 
to smaller scales. 

The anisotropic power spectrum we have constructed 
means the turbulence along different parallel rays is less 
correlated than the turbulence along two rays at the 
same separation when the power spectrum is isotropic. 
If that's the case then the transition probabilities for the 
neutrinos travelling along those two rays should also be 
less correlated and below we quantify the decrease. 



III. RESULTS 

Now that we have the random fields to model the 
turbulence we are all set to generate turbulent profiles 
and send neutrinos and antineutrinos through them. To 
achieve higher efficiency we follow six neutrinos and six 
antineutrinos simultaneously through every realization of 
the turbulence with one neutrino and one antineutrino 
emitted at a; e {0, lO'^, 10^ 10^ lO'^, 10^} cm. Each time 
we generate a new realization we end up with a different 
set of transition probabilities so by repeating the calcula- 
tion many times - in our case a minimum of one thousand 
times but often much larger - we can create an ensemble 
of transition probabilities of size N from each emission 
point. Once we have our ensemble we can then go ahead 
and compute means {Pij(x)), variances Vij{x), and, of 
course, correlations 



Pi] {Sx) 



{P^j{x)P,jix + Sx)) - {P^jjx)) {P,j{x + Sx)) 
\/Vij{x) Vij{x + Sx) 



(11) 

The correlation of the antineutrino transition probabil- 
ities will be denoted as pij. In the large N limit the 
error on the correlation is expected to be Cp = (1 ~ 
p^)/\/N — 1. Combining the results from the six emis- 
sion points we can form fifteen separations Sx so fifteen 
correlations but two points must be remembered: first, 
groups of them will cluster e.g. we will have a value for 
the correlation at Sx = 10 km but also two more at 
Sx — 9 km and Sx = 9.9 km and second, these groups of 
transition probability correlations are themselves corre- 
lated - half the data in each correlation value is the same 
for all members of the cluster. Nevertheless these clus- 
ters are useful because they serve as a consistency check 
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FIG. 4: The frequency distribution of the transition 
probability P23 for each of the neutrino emission points 
X. From bottom to top the emission points are x — 
0, lO"*, lO'^, 10**, 10'^, 10* cm. The turbulence amplitude 
is set to C* = 30%, we used Nk. = 50, Na = 5 for the 3-D tur- 
bulence field generator, and the neutrino mixing parameters 
used are those given in the text with sin^ 2^13 = 0.1. 



- we should expect the results to be similar within each 
cluster - and also they give us an indication if the error in 
the results are comparable to the expected, large-N error 
CTp given above. 

We also need to specify the neutrino mixing parame- 
ters we have used. The hierarchy will be set to normal 
and we shall comment on how our results translate to the 
inverted hierarchy. As discussed, the neutrino energy will 
be fixed aX E = 25 MeV, typical of supernova neutrino 
energies and we shall set the neutrino mixing parame- 
ters to be Sml2 = 8 x IQ-^eV^, Sm^^ = 3 x IQ-^eV^, 
sin^ 2012 = 0.83, and sin^ 26*23 = 1- The recent measure- 
ments of the last mixing angle 613 by T2K , Double 
Chooz \a§\, RENO and Daya Bay ^ are all in the 
region of 6*13 « 9°. We shall adopt this value for the ma- 
jority of this paper but this result is sufficiently new that 
we shall show on occasion results with multiple values of 
^13 in order to put this result in context. 

Finally, the turbulence amplitude will be allowed 
to vary but we shall focus upon larger values. With the 
measurement of a large value of 6*13 the turbulence effects 
are negligible for amplitudes of order ^ 1% [4l| . 




FIG. 5: The mean of the transition probability Pn (circles), 
P23 (squares) and P31 (triangles) of the neutrinos emitted at 
a; = as a function of the parameter A^^ keeping the ratio 
Nk/Nd fixed at Nk/Nd = 10. The error bars are not the 
error on the mean but rather the standard deviation of the 
samples. The turbulence amplitude is set to C* = 30% and 
sin^ 26'i3 = 0.1 



A. The point source statistics 

Before we show our results for the correlation of the 
transition probabilities as a function of the emission sep- 
aration, we consider first the statistical properties of the 
ensembles for each emission point. In addition to being 
interesting in their own right and useful as a reference, 
these calculations allow us to test that our 3D random 
field generator is working properly because the ensem- 
bles for each point of emission should be consistent and 
independent of x. 

In figure ^ we show the frequency distribution of P23 
for the six emission locations x using the mixing para- 
maters given above, = 30%, Nk = 50, Nd — 5 and 
sin^(20i3) = 0.1 and the isotropic power spectrum. The 
sample size is A^ = 3265 for each emission point. In each 
panel of the figure the reader will observe that the transi- 
tion probability is almost uniformly distributed - there is 
a slight decrease in the frequency of higher values of P23 
- but, more importantly, there is no observed trend with 
X. A closer inspection of figure ^ also hints at some 
correlation: the bottom few panels of the figure are very 
similar. We have reproduced this calculation for other 
choices of the Nk and A^^ paramaters. The results are 
shown in figure (O where we plot the mean values and 
standard deviations of Pn, P23 and P31 for ensembles of 
neutrinos emitted at a; = as a function of the param- 
eter Nk keeping the ratio Nk/Nd fixed at Nk/Nd — 10. 
There is no discernable trend with A^^ and we march on 
confident that setting A^^ = 50 and Nd = 5 does not bias 
our results. 

We now allow the values of C^, and ^13 to float and con- 
sider both the isotropic and anisotropic power spectrum. 
The evolution of the transition probability means as a 
function of C^, for the two power spectra and two choices 
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FIG. 6: Left figure, the mean of the transition probabihties P12 - top panel - P13 - center panel - and P23 - bottom panel - as a 
function of C* for neutrinos emitted from a single point. The right figure is the mean of the distributions for the antineutrino 
transition probabilities P12 - top panel - P13 - center panel - and P23 - bottom panel as a function of C* for antineutrinos 
emitted from a single point. In all panels the curves correspond to either sin^ 26'i3 = 4 x 10~* (squares) or sin^ 2&13 = 0.1 
(circles). The solid symbols denote our use of an anisotropic power spectrum, the open symbols to an isotropic power spectrum. 



of are shown in figure There are many inter- 

esting trends discussed in detail in Kneller & Mauney 
(4]| . Large amplitude turbulence works its way through 
to affect every mixing channel, not just the H resonance 
channel P23, as promised so that by = 0.5 we ob- 
serve (P12) ~ 20%, (P13) ~ 10%, (P23) 50%, and 
(P12) - 20%, (P13) ~ 5%, (P23) 1%. To put this in 
context, in the absence of turbulence all these transition 
probabilities are zero when ^13 = 9°. The only neu- 
trino mixing channel with reasonable sensitivity to ^13 
is the H resonance channel P23 and even then the dis- 
parity in (P23) at Ci, ~ 0.1 disappears by ~ 0.3. In 
contrast the antineutrinos are very sensitive to ^13 even 
at large turbulence amplitudes: the expectation value for 
P13 varies by a factor of ^ 2 when 613 is changed from 
sin^ 26'i3 = 4 X 10-4 to sin^ 26*13 = 0.1, P13 and P23 on 
the other hand change by ~ 1 — 2 orders of magnitude 
between the same limits. 

While these trends are interesting, the purpose of 
figure ([nj is to compare the use of the isotropic and 
anisotropic power spectra. Except for the H resonance 
mixing channel P23, the isotropic power spectrum gives 
values of (Pij) which are smaller than the anisotropic 
spectrum. The neutrinos are more sensitive to the tur- 
bulence when the power spectrum is anisotropic because 
the neutrinos are sensitive to the amplitude of the turbu- 



lence modes of order the neutrino oscillation wavelength 
[ssl l46j which is typically in the range of ~ 10 km in the 
H resonance region. The anisotropic spectrum removed 
all power for the fluctuations in the radial direction at the 
long wavelengths above l/fc* - which is of order 10^ km 
in our calculation - and to compensate we needed to in- 
crease the power on the smaller wavelengths which means 
and effective increase of their amplitude. In fact we al- 
ready know the exact amount the amplitude is effectively 
increased because we pointed out the 1/a factor that 
appears in the one dimensional power spectrum in the 
isotropic case compared to the one-dimensional spectrum 
derived from the anisotropic turbulence. Thankfully, our 
expectations are confirmed by figure (jH]) because the in- 
crease of all the mixing channels except P23 is on the ex- 
pected scale of a. The isotropy scale paramater fc/, which 
sets the scale in the radial direction below which the tur- 
bulence is isotropic, does not play a role for these point 
source statistics. The one-dimensional power spectrum 
along the radial direction is independent of the isotropy 
scale fc/ which can be seen when comparing equations 
(m and ((HI . So if the one-dimensional power spectrum is 
independent of kj then the effect of switching the power 
spectrum from isotropic to anisotropic is solely due to the 
removal of radial long- wavelength fluctuations. The tran- 
sition probability P23 behaves slightly differently but is 
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FIG. 7: The correlation of the transition probabihties through isotropic turbulence of various turbulence amplitudes as a 
function of the distance between emission points 5x. From top row to bottom the correlations are for P12, P13, P23, P12 and 
Pi3- The turbulence amplitudes are C* = 10% (left column), C* = 30% (center column) and C* — 50% (right column). The 
values of 613 are sin^ 2^13 = 4 x 10~* (squares joined by a solid line), sin^ 2^13 = 10^'' (triangles joined by a dot-dashed line), 
sin^ 2^13 = 4 X 10^"* (diamonds joined by a double dot-dash line) and sin"^ 2613 = 0.1 (circles joined by a dashed line). 
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entirely consistent with the understanding of the effects 
in the other channels. At smaller amplitudes and the 
smaller value of ^13 there is no effect of the power spec- 
trum switch upon (P23) because the depolarization limit 
has been reached. At the larger mixing angle depolariza- 
tion has not achieved and switching the power spectrum 
leads to the effects as seen in P12 and P13 . The two-flavor 
depolarization limit is reached for the sin^ 26*13 = 0.1 case 
when C+ ~ 30%. At around this same turbulence ampli- 
tude there begins the shift to three-flavor depolarization 
where (i-23) — 1/3. Whatever the mixing angle used, 
we see that the mean value {P23) as a function of 
using the anisotropic spectrum begins the transition at 
smaller C* than the same calculation using the isotropic 
spectrum because of the increased amplitude of the small 
scale fluctuations in the former case. 



B. The correlation through isotropic turbulence 

We now turn to the correlation of the transition proba- 
bilities as a function of the distance between the emission 
points and consider first the case of the isotropic power 
spectrum. Our result for the correlation of the transition 
probabilities, except P23, as a function of the separation 
Sx at various values of ^13 and turbulence amplitudes d, 
is shown in figure ([7]). P23 is excluded is because it is 
difficult to calculate its correlation reliably. What one 
notices immediately about the results are that pi2, P13, 
P12 and pi3 all show little sensitivity to either ^13 or 
- which is in contrast to figure The reason for the 
lack of sensitivity of these correlations to 613 and is 
explained by the exponential distributions these transi- 
tion probabilities possess. Both the turbulence ampli- 
tude and the mixing angle simply 'rescale' the ensemble 
of transition probabilities and, as equation (jlip shows, 
this rescaling cannot alter the correlation. One also sees 
that the correlation of all these transition probabilities is 
high, > 0.5, for all separations 6x < 100 km. 

In contrast the correlation of P23 is sensitive to both 
013 and C*. When is of order ~ 10% the sen- 
sitivity to 6*13 arises because the distributions of P23 at 
the different mixing angle choices are very different: for 
sin^ 2013 = 4 X lO^'' the distribution is uniform, for 
sin^ 2013 = 0.1 it is strongly skewed to small values of 
P23- As Ci, increases the sensitivity disappears because 
the distributions at each value of ^13 become similar: 
this is the same behavior seen in figure ([B]). Finally, 
for = 10% the currently preferred value of ^13 gives 
greater correlation at a given seperation than smaller val- 
ues of 013. The correlation p23 is high for Sx < 10 km, 
a scale of order the proto-neutron star diameter, at for 
= 10% and decreases rapidly as increases. For 
C-i, > 0.3 the transition probability P23 of two neutrinos 
emitted from points on the proto-neutron star separated 
by a distance greater than Sx > 1 km are essentially 
independent. 




6x [cm] 

FIG. 8: The correlation of the transition probabilities P12 
- top panel - P13 - center panel - and P23 - bottom panel - 
through anisotropic turbulence as a function of the separation 
between neutrino emission points. The turbulence amplitude 
is set at = 30% and sin^ 2613 =0.1. In each panel the cor- 
relation of the transition probabihties through the isotropic 
turbulence is shown as the solid hne. The other curves in each 
panel correspond to different values of the ratio // — ki/ki,: 
fi ~ 1 are squares joined by long dashed lines, // = 10 are 
triangles joined by dash-dot lines, and // = 100 are diamonds 
joined by dot double-dash lines. The error bars on each data 
point are estimated using the large N limit prediction. 



C. The correlation through anisotropic turbulence 

The change to the mean point source transition proba- 
bilities when switching to an anisotropic power spectrum 
is both understandable and measurable but, overall, the 
effects are small and of the order of factors of a i.e. the 
amplitude by which the small scale fluctuations in the 
anisotropic spectrum increased in amplitude compared 
to the isotropic spectrum. That insensitivity no longer 
holds when we examine the correlations of the transi- 
tion probabilities because these quantities are functions 
of the isotropy scale paramater kj. The correlations of 
the transition probabilities as a function of the separa- 
tion between the emission points is strongly sensitive to 
the amount of turbulence power in the perpendicular di- 
rections and increasing kj relative to the fixed scale fc^f 
shifts the power from long wavelength, small kx and ky, 
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IxlO'* 1x10^ 1x10^ Ixio"* 1x10^ 1x10^ Ixio"^ 1x10^ 1x10' 

5x [cm] 5x [cm] 5x [cm] 



FIG. 9: The correlation of the transition probabilities P12, P13, P23, P12 and P13 as a function of the separation between 
emission points Sx. The mixing angle 613 was set at sin'^ 2^13 = 0.1. The left column of panels is for a turbulence amplitude of 
C* = 10%, the central column for C* = 30%, and the rightmost column is C* — 50%. In each panel the squares joined by the 
solid lines are for // = 1, the triangles joined by dash-dot lines are // = 10, and the diamonds joined by dash-double dot lines 
are // = 100. 
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to much shorter wavelengths, as shown in figure ([3]). The 
effects of introducing and varying fi — kj / k^, are shown 
in figure ([8]) for the case = 0.3 and sin^ 26'i3 = 0.1. 
For fj — 1 the difference between the isotropic and 
anisotropic power spectra are minimal but at larger ra- 
tios of the two scales the correlation at some given sep- 
aration Sx drops noticeably in all three channels though 
the reduction in the correlation of P12 and P13 is not 
as severe as that for the transition probability ^23- It 
is still the case that the correlation of P12 and P13 at 
typical proto-neutron star radii of (5a; ^ 10 km is larger 
than 0.5 if //< 10 for this particular mixing angle choice 
and turbulence amplitude. Pushing even more power to 
smaller scales would lead to minimal correlation of these 
two transition probabilities. In the case of P23, the cor- 
relation at (5x ~ 10 km is already small for this mix- 
ing angle and turbulence amplitude even in the isotropic 
and fi — 1 cases so pushing more power of the fluctua- 
tions in perpendicular directions to smaller wavelengths 
completely removes the correlation of P23 over the proto- 
neutron star radial scale. 

If we now vary the turbulence amplitude we generate 
figure (jH]) which is, again, for a mixing angle of sin^ 26*13 = 
0.1. Examining the results we quickly observe the same 
general trends with changes in // as seen in figure ([5]): 
increasing // reduces the correlation with P23 affected to 
a greater degree than pi2, pis, pi2 and P13. Likewise, 
the trends seen in figure (O for changes in are also 
reproduced. 

IV. SUMMARY AND CONCLUSIONS 

Supernova turbulence and its effects upon both the 
flavor composition of neutrinos that pass through it and 
their correlations depend upon many numerous param- 
eters one needs to introduce to describe the turbulence. 
All affect the result and here we try to succinctly sum- 
marize our results. For a neutrino energy of 25 MeV and 
using a supernova density profile taken from a simulation 



4.5 s post-bounce, we find in a normal hierarchy that the 
correlation of the H resonance mixing channel transition 
probability P23 as a function of the emission separation 
Sx drops considerably as Ci, increases for both the cases 
of isotropic and anisotropic turbulence. If the turbulence 
amplitude is of order ^ 0.1 then the correlation of the 
transition probability P23 for neutrinos emitted from op- 
posite sides of the proto-neutron star, i.e. separated by 
^ 10 km, is marginal for the isotropic spectrum and for 
the anisotropic only when // < 10. For > 0.3 there 
is essentially no correlation of the H resonance transition 
probabilities for neutrinos emitted from opposite sides of 
the proto-neutron star. At these amplitudes the turbu- 
lence along parallel trajectories separated by ~ 10 km is 
just too different to permit any correlation of this tran- 
sition probability in the supernova neutrino burst signal. 
If we switch to an inverted hierarchy then it will be the 
transition probability P13 which behaves this way. 

In contrast, the correlation of the transition probabil- 
ities P12, P13, P12 and P13 in a normal hierarchy as a 
function of emission separation is largely independent 
of and 613. The correlation decreases as the ratio 
// = ki/ki, increases but remains significant for sepa- 
rations of order the proto-neutron star radius even for 
// ^ 100. When switching to an inverted hierarchy the 
mixing channels which behave this way are Pi 2, P23, Pi 2 
and P23. These mixing channels, particularly P12 and 
Pi 2, are the most promising for observing fiavor scin- 
tillation assuming the energy resolution of our neutrino 
detectors does not wash out the effect and the temporal 
correlation remains high. 
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